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ABSTRACT 



Aims. On 13 June 1998, the TRACE satellite was fortuitously well placed to observe the effects of a flare-induced EIT wave in the 
corona, and its subsequent interaction with coronal magnetic loops. In this study, we use these TRACE observations to corroborate 
previous theoretical work, which determined the response of a coronal loop to a harmonic driver in the context of ideal magnetohy- 
drodynamics, as well as estimate the magnetic field strength and the degree of longitudinal inhomogeneity. 

Methods. Loop edges are tracked, both spatially and temporally, using wavelet modulus maxima algorithms, with corresponding loop 
displacements from its quiescent state analysed by fitting scaled sinusoidal functions. The physical parameters of the coronal loop are 
subsequently determined using seismological techniques. 

Results. The studied coronal loop is found to oscillate with two distinct periods, 501 ±5 s and 274±7 s, which could be interpreted as 
belonging to the fundamental and first harmonic, or could reflect the stage of an overdriven loop. Additional scenarios for explaining 
the two periods are listed, each resulting in a different value of the magnetic field and the intrinsic and sub-resolution properties of 
the coronal loop. When assuming the periods belong to the fundamental kink mode and its first harmonic, we obtain a magnetic field 
strength inside the oscillating coronal loop as 2.0±0.7 G. In contrast, interpreting the oscillations as a combination of the loop's natural 
kink frequency and a harmonious EIT wave provides a magnetic field strength of 5.8±1.5 G. Using the ratio of the two periods, we 
find that the gravitational scale height in the loop is 73+3 Mm. 

Conclusions. We show that the observation of two distinct periods in a coronal loop does not necessarily lead to a unique conclusion. 
Multiple plausible scenarios exist, suggesting that both the derived strength of the magnetic field and the sub-resolution properties of 
the coronal loop depend entirely on which interpretation is chosen. The interpretation of the observations in terms of a combination 
of the natural kink mode of the coronal loop, driven by a harmonic EIT wave seems to result in values of the magnetic field consistent 
with previous findings. Other interpretations which are realistic, such as fundamental mode/first harmonic and the oscillations of two 
sub-resolution threads, result in magnetic field strengths that are below the average values found before. 

Key words. Magnetohydrodynamics (MHD) — Sun: corona — Sun: magnetic fields — Sun: oscillations 



1. Introduction 

Sudden energy releases in the solar atmosphere are known to 
generate large-scale global waves propagating over long dis- 
tances. The energy stored in these waves can be released by tra- 
ditional dissipative mechanisms, but might also be transferred to 
magnetic structures that come into contact with global waves. 
In this context, EIT waves generated by coronal mass ejections 
(CMEs) and/or flares could interact with coronal loops, result- 
ing in the generation of kink modes, i.e. oscillations that ex- 
perience periodic movement about the loop's symmetry axis. 
These generated loop oscillations have been used as a basic 
ingredient in one of solar physics most dynamically expand- 
ing fields, namely coronal seismology. Depending on the ori- 
entation of the studied wave, coronal seismology can be further 



[Baneriee et £ 


il. 20071 Verth et al J 20071) 


(IBallai et alJ 


2005; 


BaUail 2007b. Althou 



mology appear to be two different aspects of coronal physics, in 
reality they are closely related. 

Coronal seismology is based on theoretical relations (called 
dispersion relations) linking plasma parameters, such as the 
plasma density, to wave parameters, such as the wave frequency, 
in a precise way. In general, plasma parameters are determined 



from wave parameters, which themselves are determined obser- 
vationally. The dispersion relations for many simple (and some 
quite complicated) plasma structures under the assumptions of 
ideal magnetohydrodynamics (MHD) are well known; they were 
derived long before accurate EUV observations were avail- 
able (Edwin & Roberts 1983; Roberts et al. 1984) using simpli- 
fied models within the framework of ideal and linear MHD. 
Although the accurate interpretation of many observations can 
be difficult because of the insufficient spatial and temporal reso- 
lutions of present satellites, considerable amounts of information 
about the state of the plasma, and the structure and magnitude of 
the coronal magnetic field, can still be obtained. 

Most of the previous studies of waves and oscillations in 
coronal loops have focussed on the (local) oscillation itself, 
without considering the more global cause of the oscillation. One 
of the first studies where the nature of a local loop oscillation was 
investigated in terms of the ty pe of driver (global EIT waves) was 
the study by Ballai et alJ (l2008l) . These authors analyzed the pat- 
tern of possible oscillations (treated as an initial value problem) 
that can be recovered in a simplified loop model under the influ- 
ence of, e.g. harmonic drivers interacting with loops. Their main 
conclusion, which is used in the present study, was that for a cer- 
tain range of driver periods, the loop supports oscillations whose 
period is a superposition of the driver's period and the natural 
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period of the loop. A sim ilar analysis of a 2D loop was carried 
out bv lSelwaetall(l2010h . 

The aim of this paper is to investigate, using TRACE EUV 
observations, the characteristics of a coronal loop that oscillates 
under the influence of an external driver The special circum- 
stances of this event reside in the possibility of connecting the 
loop oscillation to a global BIT wave, presenting the phenomena 
of loop oscillations in a much wider context. 

2. Observations and data analysis 

In contrast to the SOHO/EIT instrument, the TRACE/EUV 
imager has good temporal and spatial resolution, suitable for 
observing waves in the MHD domain. However, the field-of- 
view of TRACE is merely 8.5 minutes of arc, compared to 
SOHO/EIT's full disk images of 45 minutes of arc. Thus, ob- 
taining TRACE data close to the source of an EIT wave re- 
quires a small measure of good fortune; one such instance oc- 
curred on 13 June 1998. On this day, a coronal wave was initi- 
ated just south of the field of view of TRACE at approximately 
15:23 UT by a GOES C2.9 flare. The maximum flare intensity 
occurred between the eastern part of the main negative polar- 
ity and the western part of the main positive polarity regions 
of AR 8237 (S25W04 in the solar heliographic coordinate sys- 
tem) at 15:33 UT, where a filament was also observed in Ha. 
Furthermore, this event was associated with a halo CME ob- 
served in whit e light by L ASCO C2 with a leading edge moving 
at 190 kms"' (lDelannee ''2000) . A limited study of this event has 
already been carried out by Wills- Davev & Tho mpson (1999'), 
with the authors interpreting the disturbance as a fast magnetoa- 
coustic wave. 

TRACE images in the 171 A and 195 A bandpasses ac- 
quired during 15:00 - 19:00 UT show the progress of the wave 
(activity over an extended area of the field-of-view), and from 
the SOHO/EIT instrument a con comitant EIT wave is visible 
(IWiUs-Davev & Thompsonlll999 ). The TRACE data also reveal 
how the EIT wave front causes the positioning of some coronal 
loops to be disturbed. The particular loop under investigation is 
situated in the south-eastern part of Fig. [T] with a zoom into the 
region of interest in Fig.|2] 

Using the standard calibration/correction packages available 
in the SolarSoft library, the TRACE data was processed in the 
usual way: spikes from cosmic ray hits were removed, data 
frames were normalized to a constant exposure time, and the 
image sequence was de-rotated. The TRACE satellite entered 
the Southern Atlantic Anomaly at approximately 15:50 UT, and 
images up to this point clearly display a propagating wave front. 

2. 1 . Temporal analysis of loop edge 

Th e characteristics of a global EIT wave were studied in detail 
bv lBallai et al.l (l2005h . who demonstrated that the global distur- 
bance is a wave with an average periodicity of approximately 
400 s. To accurately probe the eff'ects of this wave front on more 
rigidly fixed coronal loop structures, the loop edge was tracked 
spatially and temporally using a wavelet modulus maxima tech- 
nique described by iJess et al ] (120081) . Since the 195 A response 
function is somewhat lower than that of the 171 A passband, 
longer exposure times are required to provide the same signal as 
the 171 A channel (a 41 s average exposure time at 195 A, ver- 
sus 24 s at 171 A images). As a result, rapidly moving structures 
may be blurred in the 195 A passband because of their substan- 
tially longer exposure times. Furthermore, the 171 A images are 
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Fig. 1. TRACE field-of-view at 171 A immediately prior to the 
arrival of the EIT wave front. The white box indicates the coronal 
loop under investigation in the present study. Axes are in solar 
heliocentric arcseconds. 



at full resolution (1 pixel = 0.499"), whereas the 195 A data are 
2x2 binned, resulting in a pixel size of 0.998". Through use of 
blin d iterative d econvolution techniques, Golub et al. ( 1999) and 
ILin et al.l (l200lb determined the absolute spatial resolution of the 
TRACE instrument to be 1.25". Since the 171 A data is not de- 
graded by spatial binning, it enables us to monitor small-scale 
disturbances more clearly in coronal loop structures. As a re- 
sult, only images taken in the 171 A bandpass could be tracked 
with an average cadence of 90 s. First, a segment of the loop 
edge, between the loop footpoint and apex, was passed through 
a Laplacian filter to sharpen the images and make the edge of 
the structure more pronounced. Next, a lower intensity thresh- 
old was applied to the data to remove contributions from the 
underlying quiet Sun. A threshold of the background mean flux 
plus 4.5 cr was used, which produces the image displayed in the 
upper-left panel of Fig. |2] To emphasize feature edges and re- 
move shallow intensity gradients, we used a binary format for 
the feature mapping. All pixels of values above the lower flux 
intensity threshold defined above were assigned a value of '1'. 
Those pixels which lie below the threshold are assigned a value 
of '0', producing the binary map shown in the upper-right panel 
of Fig. |2] To assist the analysis, each image was rotated by 45° 
clockwise until the southern edge of the loop structure was par- 
allel with the X-axis. This caused loop displacements to be con- 
fined to the direction of the y-axis, as a function of the loop 
length (x-axis). TRACE images acquired between the initiation 
of oscillatory behaviour at 15:29 UT, and 15:50 UT when the 
instrument entered the Southern Atlantic Anomoly, were binary 
mapped for e dge detec tion and tracking using the methodology 
described by Jess et al.1 (l2008h . 

Utilizing the resulting loop-edge positions, the displacement 
of the loop segment as a function of time from its quiescent 
state was established for 30 pixels (wl 1000 km; see dashed lines 
in Fig. |2]l along the loop edge between the footpoint and apex. 
Owing to the minimal shifts in the loop equilibrium position with 
time, and the need to minimise the number of free fitting param- 
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Fig. 2. Zoomed images of the TRACE field-of-view contained within the white box in Fig. [T] The upper-left image displays a 
true inten sity image, while the upper-right panel reveals the corresponding loop-edge binary map used in the tracking algorithm of 
iJess et all (12008). Locations of the loop footpoints are highlighted with green boxes in the left panel, giving an approximate loop 
length, L — 72.9±2.8 Mm. An 11000km segment, incorporating a brightly emitting section of the loop edge, which is used in the 
search for higher harmonics is contained within the green dashed lines plotted in the upper panels. The lower-left panel displays 
an emission measure map, derived from near-simultaneous 171 A and 195 A TRACE images, with the corresponding color scale in 
units of cm"^. Black boxes and dashed lines represent characteristics identical to those described above for the upper-left panel. The 
loop arcade intensity (+), averaged between the green dashed lines in the upper-left panel, is plotted in the lower-right panel as a 
percentage above the quiescent background. A Gaussian fit to the coronal loop under investigation is displayed using a solid black 
line. Green shading represents loop intensities greater than the fitted half-maximum, while cross-hatched red shading indicates 
regions of the loop structure with intensities below the half-maximum. Arrows mark the full-width at half-maximum, chosen to 
represent an estimate of the loop depth, found here to be =s8 pixels, or siSOOOkm. 



eters, no long-term trends were subtracted from the original time 
series. The loop-edge displacement curve was then fitted using 
sinusoids of varying frequency and amplitude. The best fit si- 
nusoid was then subtracted from the original data to produce a 
residual time series. The residual data, again, was modelled with 
the best-fit sinusoid function. To determine the success of the fit- 
ting process, Kolmogorov-Smirnov and cross-correlation statis- 
tics (iPress et al...l992; .Christian et al... 1998: .Jess et al...2007a.bi) 



were determined for each stage of the sinusoid fitting process, 
as well as for the combined fit with respect to the original 
time series. An example of this process is displayed in Fig. [3] 
where Kolmogorov-Smirnov statistics for the 501 s, 274 s, and 
combined fits are 88%, 91%, and 92%, respectively. Cross- 
correlation statistics for the same fitting functions are 90%, 74%, 
and 98%, respectively, indicating that these sinusoid functions 
recreate the observed data with a high degree of precision. 
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Fig. 3. Top Left: A typical transverse displacement of one pixel along the loop segment (+), joined together using a spline fit to 
assist visualisation of the motion path (solid line). Top Right: The original transverse displacement (+), overplotted with a 501 s 
period sinusoidal signal (dashed line). Bottom Left: The residual (501 s period sinusoid subtracted from the original motion path) 
transverse displacement (+), overplotted with a 274 s period sinusoidal signal (dashed line). Bottom Right: The original data points 
(+) and spline fit (solid line) shown in the upper-left panel, overplotted with the combined 501/274 s period sinusoidal fits. These 
combined fits, with respect to the original time series, result in Kolmogorov-Smknov and cross-correlation statistics of 92% and 
98%, respectively, indicating a high degree of accuracy. 



To validate this result, the process of loop isolation, tracking, 
and sinusoid-fitting was repeated for image sequences that had 
been spatially degraded using 2x2, 3x3, 5x5, and 7x7 binning 
intervals. While spatial binning has no effect on a robust peri- 
odic signature, it will help us to reduce any artefacts present in 
the data, such as those caused by detector noise and/or cosmic 
ray hits. These independent analyses yielded consistent results, 
for which periodicities overlapped with those established above, 
implying that the oscillation occurs over a large spatial scale, 
an d is unaffected by varying smo othing intervals. As described 
in IVan Doorsselaere et alJ (l2007l) . repeating this procedure for, 
and averaging over, many loop positions and binning intervals 
(30 loop positions x 5 binning values =150 independent mea- 
surements in this instance) allows the associated errors to be 
drastically reduced. Following a similar methodology, we derive 
best- fit periodicities of 501+5 s and 274+7 s, both of which are 
above the limiting Nyquist period (»;180 s) of the dataset. 

2.2. Properties of the coronal loop 

The fact that two different periods have been recovered from 
the same structure allows us to uncover properties of the loop 
that cannot be directly measured. We use a simple semicircular 
model for the loop, where the loop length, L, is given by 7tD/2, 
where D is the distance between footpoints. As the TRACE field 
of view is nearly at the solar disk centre, and the loop does not 
appear to be tilted at a significant angle, this model is likely to be 



accurate. The distance D is measured to be 128+5 TRACE pix- 
els, corresponding to an absolute diameter of D - 46.4+1 .8 Mm. 
Thus, we are dealing with a relatively short loop, with a corre- 
sponding estimated loop length of L = 72.9+2.8 Mm. 

The electron number density of the loop (and its environ- 
ment) can be estimated using a differential emission measure 
technique by comparing the intensities of the loop detected at 
different wavelengths. This filter-ratio method is based on the 
assumption that a structure seen in both 171 A and 195 A band- 
passes will be isothermal, with an associated temperature falling 
between the two peak temperatures of the TRACE filter response 
functions (Aschwanden et al. 2 000i) . Using near-simultaneous 
171 A and 195 A TRACE ima ges, an d following the methodol- 
ogy described in i AschwandenI (120041) . we calculate the emission 
measure, EM, as EM = (4.0+0.2) x 10^^ cm"^, averaged over 
the 30-pixel loop segment used in Section ITTI (lower-left panel 
of Fig. |2]i. This value represents the emission integrated over the 
entire line-of-sight. The relationship between emission measure 
and electron number density, n^, in a source plasma structure of 
depth, d, is given by 

j EM(x,y) 

ne(x,y) = ^J , (1) 

where x and y are the spatial coordinates (lAschwandenll2004h . 
Assuming a cylindrical loop structure, and using the full-width 
at half maximum to estimate the loop cross-section, we estimate 



Istvan Ballai et al.: TRACE observations of driven loop oscillations 



5 



the loop depth to be 8+1 pixels, or (29.2+3.6) x 10^ cm (lower- 
right panel of Fig. This allows us to calculate the electron 
number density in the loop as n^. = (3.7±2.3)x 10*^ cm"-'. Using a 
quiescent plasma electron number density (the region outside the 
loop) of 1 + 0.5 X 10^ cm"^* in line with previous estimates (see, 
e.g. Fludra et al. 1999, Ugarte-Urra et al. 2002), we can deter- 
mine the loop filling factor to be 3.7+0.4. For brightly emitting 
coronal loops, it can be assumed that all of the emission measure 
comes directly from the loop itself. However, for relatively dim 
loops, similar to the one under investigation here, a significant 
portion of the emission will come from plasma in the immedi- 
ate vicinity of the loop. As a result, the derived electron number 
densities will represent the upper boundaries of actual values. 

2.3. Application of coronal seismology 

Coronal seismology, initiated after high-resolution observations 
of waves and oscillations in coronal structures became available, 
can provide relatively accurate estimates of quantities that can- 
not be measured directly or are below the instrument's resolution 
limit. Using a combination of observations and theoretical mod- 
els, answers to the problems of the magnitudes and structuring 
of coronal magnetic fields, density scale-heights, etc ., have since 
been recovered d Nakariakov et al. 1999; Andries et al. 1 120051 
l2009tlVerth et aLlEoOStlRuderman & Erdelvi 2009). 

In this context, one of the most promising methods for in- 
vestigating the internal structure of coronal loops is based on 
the P\/P2 ratio of the fundamental (Pi) to first harmonic (P2) 
modes of an oscillation. In a homogeneous, non-stratified at- 
mosphere the period ratio is 2. However, longitudinal varia- 
tions in the properties of the plasma (i.e. variations along the 
loop) can make this ratio to deviate either below or above 2. It 
is believed that magnetic effects produce a ratio that is greater 
than 2 because of the magnetic dive rgence caused by changes in 
the cross-sectional area of the loop (iDe Moortel & B radv 2007; 
lO'Shea et al.ll200'7h . In this case, the magnetic efl^ect on the pe- 
riod of oscillation dominates any other eff'ects, such as den- 
sity stratification. In contrast, if the ratio of the two periods 
is less than 2, then the change is attributed to density strati- 
fication effects along the coronal loop. Several examples have 
been published where the P1/P2 ratio is within the interval 
1.5 - 1.9 (McEwan et al. 2006; Van Doorsselaere et al. 200 
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Ruderman et al. 2008; Verth & Erdelvill2008t iVerth et al.ll2008l 
Andries et al. (2005) show how a one-to-one relationship exists 
between the P1/P2 ratio and the density scale-height. The de- 
pendence of the period ratio on L/nH, where H is the density 
scale-height, is sh own in Fig.Bl 

According to Andries et al.l (2005), the dependence of the 
period ratio on the density stratification does not depend on 
other loop parameters or the magnetic field strength. As a re- 
sult, the same dependence can be demonstrated for a vari- 
ety of density profiles generated inside the same coronal loop. 
For the coronal loop under current investigation, we find that 
P\/P2 = 1.82+0.02, corresponding to a density scale height 
H - 73±3 Mm. This value i s consistent with those found by 
Doorsselaere et a n (l2007h . 

Assuming that the loop is in hydrostatic equilibrium, the den- 
sity scale height is directly proportional to the temperature, T, 
according to the equation 



H = 47 ( ) ,(Mm). 

\mkI 



(2) 




Fig. 4. The period ratio, P1/P2, corresponding to the fundamen- 
tal to first harmonic kink mode frequencies, plotted as a function 
of loop length, L, and density scale height, H. A P1/P2 ratio 
equal to 1.82, combined with a loop length La;73 Mm, provides 
a density scale height of H^12 Mm. 



sponse for the 171 A channel, and towards the central portion of 
the 1 9 5 A channel's temp erature response curve (lAschwandenl 
120041: iPhillips et all l2005h . Nevertheless, the visibility of the 
coronal loop in both the 1 7 1 A and 1 95 A bandpasses, in addition 
to a derived loop temperature within the TRACE instrument's 
temperature response function, suggests that the loop under in- 
vestigation is in hydrostatic equilibrium dAndries et alj|2005l) . 

On the basis of the peri od of the fundamental m ode (501 s), 
following the methods of Verwichte et al ] (I2004h . It IS possi- 
ble to calculate the kink speed, ck, of the wave to be ck - 
2L/Pfu,idamentai = 291+8.2 kms"^ The kink speed can then be 
related to the magnetic field strength, B, by the equation 



B^CK 



l^iPi + Pe) 



1/2 



(3) 



Thus, for our present analysis we obtain a loop temperature of 
T - 1.5+0.6 MK, which is within TRACE'S temperature re- 



where p, and pe are the internal and external loop densities, re- 
spectively, yu is the permeability of free space, and the magnetic 
field strength inside and outside the loop is assumed to be iden- 
tical. Thus, we estimate the magnetic field strength inside the 
coronal loop to be 2.0+0.7 G, which is towards the lower-en d of 
statistical studies compile d bv lNakariakov & OfmanI (12001 h and 
lAschwanden et al.1 (l2002h . 



2.4. Comparison with other theories 

The existence of two distinct periods in the loop under inves- 
tigation has been explained and investigated assuming that the 
two periods belong to the fundamental and first harmonic kink 
oscillations of the loop. However, other viable explanations are 
also possible. For our discussion, we first recall the results of 
IBallai et alJ (l2008h . Using a simple coronal loop model, these 
authors investigated the interaction between an EIT wave and 
a coronal loop assuming an equilibrium of forces. Their results 
show that the periods recovered in a coronal loop are always a 
combination of the period of the driver (here EIT waves) and the 
natural period of the loop. The same authors showed that the dis- 
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placement of a coronal loop (here denoted by Q(z, t)) under the 
influence of a driver, F(z, t), is given by 



5^ 



(4) 



where, as before, ck is the kink speed (propagation speed of dis- 
turbances in the loop) and is the cut-off frequency of kink- 
mode oscillations. These quantities are simply given by 



Ck 



i 



Pi + Pe 



g(Pi-Pe) 
2cK(Pi + Pe) 



where v^, and v^e are the Alfven speeds inside and outside the 
loop. Equation |4] is an inhomogeneous Klein-Gordon equation 
that was solved assuming that the foot-points of the loop are 
fixed (line-tying condition) and initially at rest. 

The event occurring on 13 June 1998 was also studied in 
detail by Wil ls-Davev & Thompsonl ( ll999h . who showed that the 
blast wave intersects the loop of interest nearly perpendicularly, 
meaning we can model the EIT wave as a harmonic driver whose 
temporal and spatial dependence is given bv lBallaietal.l(i2008h 
to be 



F(z, t)ocK [6(z - zo) -6{z-L + z,,)] e" 



(5) 



where 5{z) is the Dirac-delta function, and the constant K de- 
pends on the density of the plasma and energy of the EIT wave, 
both considered to be constant. In Equation |5j loeit is the fre- 
quency of the EIT wave, which is assumed to be a harmonic 
signal. In this study, we are concerned with temporal changes 
only, which is why the exact form of the driver is not explicitly 
given. The form of the driver given by Equation |5] reveals that 
the interaction between the loop and the EIT wave occurs in two 
points (at z - Zd and z - L - zo) simultaneously, and is placed 
symmetrically with respect to the end-points of the loop at z = 
and z-L. 

Equat ion [H wher e F(z, f) is defined by Equation |5] was 
solved bv lBallai et al.l (12008 ), with the temporal dependence of 
the displacement, Q{z, t), found to be 



sm 



loeit 



\t sm I — 



(6) 



where o)^ = a>^ + rP'n^c^^/L^ {n > 1) are the natural frequencies 

of the loop. 

According to the results of iBallai et al.l (l2008l) (their Figure 
3), the dominant (and strongest) signal in the loop is produced 
by the driver Bearing in mind the previously presented model, 
we identify the two measured periods as the two temporal de- 
pendencies of the solution given by Equation|6] Since it is visu- 
ally obvious that the loop oscillates at its fundamental frequency, 
(x)„ in Equation|6l it is replaced by lo\ corresponding to the fre- 
quency of the fundamental mode. That means that we need to 
solve the system 



An 



4n 



a>i + ojeit — T;r^ a>\ — coeit — — , 

12 1 1 



(7) 



where T\ - 501 s and T2 - 274 s. It is easy to show that the 
period of the fundamental mode is 177.1 + 3.5 s and the period 
of the driving EIT wave is 604.7 + 26.8 s, much higher than the 
period derived by Ballai et al. (2005). However, we should keep 
in mind that the result derived by Ballai et al. (2005) refers to an 
average period derived from within a large field-of-view. 



Once the period of the fundamental mode is derived, it is 
straightforward to estimate the magnetic field strength inside 
the loop. W e, again, need to use the mathematical finding of 
'Bal lai et aP (|2008) to derive the magnetic field. According to 
their analysis, the frequency of the fundamental mode is given 
by 



-4 



Lot + n^c 



(8) 



Using the previously determined values for densities and length 
of the loop, together with the derived frequency of the funda- 
mental mode, we obtain the fundamental kink speed of the loop 
to be 823 + 15kms Using the standard definition of the kink 
speed, and assuming the same magnetic field strength inside and 
outside the loop, we can obtain the magnetic field in the loop to 
be 



B^CK 



t^iPi + Pe) 



5.8 + 1.8 G 



(9) 



The values of the magnetic field strength determined using two 
different scenarios are rather different, meaning that the deter- 
mination of this fundamental quantity depends on the applied 
theoretical model used to explain the observed periods. 

As we specified earlier the interpretation of two periods ob- 
served in a loop is not unique because different scenarios can 
occur Unfortunately, given the present spatial resolution con- 
straints, it is impossible to distinguish between these models, 
meaning that all findings derived from the observations of loop 
oscillations should be treated with care. 

Another possible explanation of the two observed periods is 
connected to the limits of the spatial resolution of the TRACE in- 
strument. The two observed periods could belong to two neigh- 
boring thin loops that cannot be resolved by our observations 
moving in phase. This possibility was studied using numerical 
simulations by Lunaetal. (2008), and using theoretical mod- 
els by Van Doorsselaere et al. (2008), Roberts on et all (l2010h . 
and Robertson & RudermaaC201 1) . Assuming that the loops are 
identical, we can easily find that the ratio of d/R, where d is the 
distance between the longitudinal axes of the two loops and R 
is their radius, is 2.00 + 0.11. Thus, the physical distance be- 
tween the two loops, in units o f their radius, is 0.031 ± 0.001 . 
Furthermore, using the results o Doorsselaere et al. I (l2008h . 
we can estimate the magnetic field inside the loop to be 2.6 + 
0.4 G, a value remarkably consistent with that obtained when in- 
terpreting the two periods as belonging to the fundamental and 
first harmonic kink mode of a single coronal loop. 

Finally, the fourth possible explanation is related to the in- 
trinsic properties of the coronal loop where the second harmonic 
is not a proper eigenmode but just reflects a slightly more trian- 
gular amplitude profile than sinusoidal. In other words, the sec- 
ond period is observed as some nonlinear amplitude oscillation 
caused by a gradient in either density or magnetic field across 
the loop motion, or a coupling with neighbouring loops. 

3. Conclusions 

To carry out coronal seismology it is imperative to know the re- 
lationship between the composition of a plasma structure and the 
oscillations supported by the coronal loop. High resolution ob- 
servations enable us to accurately measure not only the magnetic 
field strength, but also the thermodynamical state of the plasma. 

The event studied here occurred on 13 June 1998, and was 
recorded by the TRACE instrument in both 171 A and 195 A 
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bandpass filters. A sudden energy release outside the field-of- 
view generated a global wave that interacted with local coro- 
nal loops. The resulting loop oscillations were identified as kink 
modes, i.e. oscillations that propagate along the coronal loop 
leading to the displacement the longitudinal symmetry axis of 
the tube. Rigorous analysis of the oscillation revealed the exis- 
tence of two distinct periods contained within a relatively short 
loop. Periodicities of 501+5 s and 274+7 s allowed us to draw 
conclusions about the magnetic field strength and the longitu- 
dinal density structure inside the loop. The magnetic field es- 
tablished here has a value towards the lower end of recent sta- 
tistical studies. However, a 30% increase in the magnetic field 
strength can be achieved if the loop filling factor is doubled. A 
major conclusion of our study is that two observed periodicities 
in a loop structure might not be a unique interpretation, often 
with derived system parameters differing substantially owing to 
the independent seismological methods. We have compared the 
observational evidence with two different theoretical interpreta- 
tions of oscillations. First, we considered that the two periods 
are the periods corresponding to the fundamental and first har- 
monic. The P\/P2 - 1.82+0.02 ratio of the two periods has 
allowed us to determine values of 73+3 Mm and 1.5+0.6 MK 
for the density scale height and loop temperature, respectively. 
Furthermore, an estimate of the magnetic field strength inside 
the coronal loop was determined to be 2.0+0.7 G. All values 
determined are consistent with previous studies involving coro- 
nal seismology, even though the coronal loop structure presented 
here is relatively confined and dim. 

A second possible explanatio n of two periods might be that 
developed by Bal lai et al] (|2008'), which states that the two pe- 
riods might belong to the driver (here an incident FIT wave as- 
sumed to be harmonic) and the visi bly clear fundamen tal mode 
kink oscillation. Using the results of i Ballai et al.l (l2008h . we ob- 
tained the period of the independent driver and kink oscillation, 
with the latter quantity allowing us to again estimate the mag- 
netic fie ld strength inside the loop structure. The theory devel- 
oped bv lBallai et al.l (l2008l) relies heavily on the FIT wave being 
described using wave theory, as opposed to a propagating mag- 
netic feature. Our observations provide additional support that 
wave theory can describe the mechanism behind FIT waves. As 
such, we feel the use of FIT waves for coronal seismology is an 
accurate and reliable approach. The magnetic field determined 
in light of this scenario is almost triple that obtained utilising the 
first method. On the basis of the magnetic field strength alone, 
and its consistency with previous estimated values, we deem that 
the interpretation of the two periods as the combination of the 
loop's natural fundamental period and the period of the harmonic 
driver, seems more plausible. Furthermore, when the two periods 
were assumed to belong to two unresolved coronal threads we 
were able to estimate both the magnetic field inside the threads, 
in addition to the distance between them. Finally, another possi- 
ble explanation of the two periods was linked to the significant 
density or magnetic field gradient along the motion path of the 
oscillating loop, or an aspect of coupling with neighbouring loop 
structures. 

Our study reveals that the periods of oscillations of a coro- 
nal loop should be interpreted with great care as a multitude of 
scenarios are plausible that would each result in different values 
of, e.g. magnetic field. Further higher resolution observational 
evidence will help us to apply the most likely interpretation of 
dynamical processes. 
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